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Abstract. 

It is well known that a particle in a periodic potential with an additional constant 
force performs Bloch oscillations. Modulating every second period of the potential, 
the original Bloch band splits into two subbands. The dynamics of quantum particles 
shows a coherent superposition of Bloch oscillations and Zener tunneling between the 
subbands, a Bloch-Zener oscillation. Such a system is modelled by a tight-binding 
Hamiltonian, a system of two minibands with an easily controllable gap. The dynamics 
of the system is investigated by using an algebraic ansatz leading to a differential 
equation of Whittaker-Hill type. It is shown that the parameters of the system can 
be tuned to generate a periodic reconstruction of the wave packet and thus of the 
occupation probability. As an application, the construction of a matter wave beam 
splitter and a Mach-Zehnder interferometer is briefly discussed. 



PACS numbers: 03.65.-w, 03.75.Be, 03.75.Lm 
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1. Introduction 

Since tlie beginning of quantum mechanics, the dynamics of a particle in a periodic 
potential is a subject of fundamental interest. The eigenenergies of such a system 
form the famous Bloch bands, where the eigenstates are delocalized states and lead to 
transport [1]. If an additional static field F is introduced, the eigenstates become 
localized. Counterintuitively, transport is dramatically reduced and an oscillatory 
motion is found, the famous Bloch oscillation. During the last decade these Bloch 
oscillations have been experimentally observed for various systems: electrons in 
semiconductor superlattices [2], cold atoms in optical lattices [3], light pulses in photonic 
crystals [4,5] and even mechanical systems may Bloch-oscillate [6]. 

If the external field is strong enough, tunneling between Bloch bands becomes 
possible, as already discussed by Zener long ago [7]. In general, successive Zener 
tunneling to even higher bands will lead to decay as the bandgaps usually decrease with 
increasing energy. This decay can be observed as pulsed output for a Bose-Einstein 
condensate in an optical lattice [8], for electrons in semiconductor superlattices [9] or 
for light in waveguide arrays [10]. In this case a description in terms of Wannier-Stark 
resonances is more suited than the common Bloch band picture [11]. 

In order to study and exploit the interplay between Bloch oscillations and Zener 
tunneling, a system of at least two bands is needed, which are well separated from all 
other bands. This is achieved by introducing another weak potential with a doubled 
period length. Due to this pertubation, the ground band splits into two minibands and 
Zener tunneling between the minibands must be taken into account. 

In this paper we discuss the dynamics in such a period-doubled potential within 
the tight-binding approximation. We assume that the ground state energies of the 
potential wells are alternately increased or decreased by the additional potential. The 
Hamiltonian then reads 

. oo r oo OO 

H = -- {\n){n + l\ + \n){n-l\) + - {~mn){n\ + Fd ^ n\n){n\, (1) 

n=— OO n=— oo n=— oo 

where \n) is the Wannier state localized in the nth potential well. Furthermore, A 
denotes the width of the unperturbed band {6 = 0), 6 gives the gap between the two 
minibands at the edge of the Brillouin zone and F is the external force. For 6 = 
this Hamiltonian reduces to the celebrated single-band tight-binding model, which has 
been studied in great detail (see [12] for a recent review). The period-doubled tight- 
binding model was introduced for the field-free case F = already in 1988 [13]. First 
pertubation theoretical results with an external field can be found in [14,15]. Numerical 
results for an analogous two-miniband system in two spatial dimensions were reported 
in [16,17]. 

This paper is organized as follows: First of all we repeat some results for the 
field free system F = 0. The spectral properties for the case F ^ are discussed 
in section IHl and the dynamics is analyzed in detail in section |3] We present an 
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algebraic ansatz for the time evolution operator and analyze the time evolution by 
general arguments. The interesting case of periodic, i.e. reconstructing, Bloch-Zener 
oscillations and thus reconstructing occupation probabilities is analyzed in detail. These 
theoretical considerations are illustrated by some numerical results. As an application of 
Bloch-Zener oscillations, a beam splitting mechanism for matter waves is described and 
simulated numerically. The possibility of constructing a Mach-Zehnder interferometer 
based on Bloch-Zener oscillations is briefly discussed. 



2. Bloch bands in the field-free case 



First of all we discuss the spectrum of the period-doubled Hamiltonian (0) for the field- 
free case F = 0, i.e. Bloch bands and Bloch waves. A straightforward calculation yields 
the dispersion relation 



^/3,K = |(-If+V5' + A'C0S2M, 



(2) 



with the miniband index /5 = 0, 1 and 7 = sgn(5), which is illustrated in figure ^ For 
6^0, the Bloch band splits into two minibands with band gap 6. 



LU 




-25 
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Figure 1. Dispersion relation of the double-periodic system in the reduced Brillouin 
zone for A = 80. The Bloch bands for 6 = (dashed curve) split into two minibands 
with gap 6 (solid lines for 6 = 18). 

The Bloch waves \xi3,k) of both bands are given by 

^00 _ 00 

^ n=— 00 ^ n=— 00 

^00 _ CO 



1X1 



n=—oo * n=— 00 



with normalization constant A^^ = vr('u^ + v^)/d and coefficents 

= A cos(/t(i) and = S + ■y ^/6^~+~A^~co^iJ^ . 



(5) 
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In the Bloch basis, the time independent Schrodinger equation with Hamiltonian ((H) 
reads 

E(xo,k|^) = ^o,«(xo,«|^) + iF^(xo,K|*) - Fd{xo,n\'i^) - M,(xi,.|*) (6) 

E{xi,.\^) = + iF^(xi,.|vl>) + M«(xo,.|^) , (7) 

where the couphng of the two bands is given by the transition matrix element 

{Xi,.'\Hzb\Xo,.) = MJ^{k - k') (8) 
with the reduced matrix element 

(1:^.) V. - {i:v.) _ iFdA6 smjKd) e^^'^ 

ul + vl 2(52 + 2A2cos2(Kd) ^ ' 

(cf. [14]). The modulus of is shown in figure |21 




Figure 2. Absolute value of the reduced transition matrix element |Mk| for A = 80 
and dF = 1. The range |Mk| > 0.5 has been cut off to increase the visibility. 

Obviously, the band transitions mainly take place at the edge of the Brillouin zone and, 
because of the delta- function in (jHJ, only direct interband transitions are possible. 

3. Spectral properties 

We now discuss the spectrum of the two-miniband Hamiltonian (Q) with an external 
field. The spectrum of the single-band Wannier-Stark Hamiltonian {6 = 0) consists of a 
ladder of equally spaced eigenvalues - the Wannier-Stark ladder. Note that beyond the 
single-band tight-binding approximation these eigenstates become resonances embedded 
into a continuum (see [11] for a review). In the case of two minibands, 5 7^ 0, the 
spectrum consists of two ladders with an offset in between, as will be shown below. 
Furthermore the dependence of the offset on the system parameters will be discussed in 
detail. Certain aspects of the spectrum were previously discussed in [14]. 
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To keep the notation simple, we introduce the translation operator 

oo 

fm= ^ \n- m){n\ (10) 

?i=— oo 

and an operator G that causes the inversion of the sign of 6 in all following terms. They 
fulfil the following commutation relations 



6 



^ E (-ir[i-(-ini^-^)(^i 

n=—oo 

oo 

+ Fd m\n — m){n\, 



'111 



S J2 and 



An eigenvector oi H with the eigenvalue E{5, A, Fd), 
H\^) = E{6,A,Fd)\^) , 



(12) 



(13) 



satisfies 



I: 



H^T2i+iG\^) 



-- {E{5,A,Fd)-2ldF}{T2i\<il)j , (14) 
= {Ei-6,A,Fd)-{2l + l)Fd}{f2i+iG\^)] . (15) 

Thus the application of the operators T21 resp. T2«+iG' to the eigenvector yields two 
ladders of eigenvectors 

H\<ila,n) = ^a,n|^a,n) with « = 0, 1 (16) 

defined by 

l^o.n) = r_2„|^) and =Gf_(2„+i)|^) (17) 

with equidistant eigenenergies 

Eo,n = E{5,A,Fd) + 2ndF , (18) 

Ei^n = E{-5,A,Fd) + {2n + l)dF . (19) 

Thus the spectrum of the Hamiltonian (P) consists of two energy ladders with equal 
spacings and an offset given by E{S, A, Fd). Furthermore it can be shown that there 
are not more than two different energy ladders (see appendix). The eigenvectors |^I/a,n) 
are displaced by 2n lattice periods with respect to I^E^a.o), where a labels the two ladders. 

The dependence of E{6,A, Fd) on the system parameters 6, A and Fd shows several 
interesting symmetries. In the case 5 = it is well known that for all values of A 
and Fd we can choose E{0, A, Fd) = [12]. In the case S ^ it can be shown that 
E{-6, A, Fd) = -E{S, A, Fd) holds, as well as E{6, - A, Fd) = E{6, A, Fd). The proofs 
of these relations are quite similar and only one of them is given here. We consider the 
operator X defined by 



X= J2 \-n){n\i-l)r 



(20) 
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which fulfils the relation 

XH{6, A, Fd) = -H{-6, A, Fd)X. (21) 
Applying —X to equation (fTI)jl yields 

-E,,„(5,A,Frf)|x|^„,„)} = if(-5,A,Frf)|x|^„,„)} 

= : E^,,n'{-S,A,Fd) . (22) 

This means that the whole spectrum is antisymmetric by means of 

- E„,„(5, A, Fd) = E^,^n'{-5, A, Fd) . (23) 

To simplify the following argumentation, we consider without loss of generality the whole 
spectrum S modulo 2dF. Thus we get 

S = {E{5, A, Fd), E{-6, A, Fd) + dF^^^^dF ■ (24) 
For —6 the spectrum therefore reads 

S = {E{-6, A, Fd), E{6, A, Fd) + dF}^,^2dF • (25) 
By equation this must be equal to 

S = {-E{6, A, Fd), -E{-6, A, Fd) - dF^^^^dr 

= {-E{6, A, Fd), -E{-6, A, Fd) + dF]^,^2dF ■ (26) 
From the equivalence of the equations and it follows that 

E{-6,A,Fd) = -E{6,A,Fd). (27) 
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Numerical results for the energy offset E{S,A,1) and the energy ladders 
£'cf^„(5, A, 1) are shown in figure El and figured In all examples we consider only the 
case Fd = 1, as all other cases can be reduced to it by a simple scaling 

A, Fd) = dFE^^^ A, 1) , (28) 

which directly follows from a rescaling of the Hamiltonian H by Fd. 

In certain limits, one can easily derive some analytic approximations for the ladder 
offset E{6,A,Fd): 

(i) For 6 ^ A the two minibands are well separated energetically, so that their coupling 
can be neglected. Hence one sets ^ in the equations © and ((7j). The energy 
shift is then given by 



/ ^/6^ + A^cos^ydy (30) 
Jo 



EiS, A, Fd) = ^V6^T^I (f , (29) 

TT 

(compare [14]) with 7 = sgn((5). Here, / denotes the incomplete elliptic integral of 
the second kind [18]. This approximation is compared to the numerical results in 
figure m 

(ii) For small 6 perturbation theory can be done, which yields 

E{6,A,Fd) = po{f-,) , (31) 
where Jq denotes the ordinary Bessel function. 

The zeros of E{6, 0,1) are 6 = 0, ±2, ±4, . . . (see figure |2I), whereas the zeros of 
E{6,A,1) for very small 6 are given by the zeros of Jo(A) (see equation ()31|)) which 
are A ^ rnc for large A. These relations indicate that the characteristic structure of 
the energy shift looks similar for large A and large 6. In figure |3] one clearly sees the 
asymmetry of E{6, A, 1) with respect to 6. Avoided crossings appear in the typical case 
of A ^ 0. 

4. Dynamics 

4.1. An algebraic ansatz for the time evolution operator 

In this section we derive some general results for the time evolution operator U. The 
single-band system 6 = can be treated conveniently within a Lie-algebraic approach 
[19]. To simplify notation, we rewrite the Hamiltonian using the operators 

00 00 
N=J2n\n){n\, L = (-l)-|n) (n| = (-1)^ , (32) 

71=— 00 n=— 00 

00 00 

K = \n-l){n\ and = \n + l){n\. (33) 
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Figure 4. Spectrum Ea,n{S) of the Hamiltonian for Fd = 1 (solid lines). The left- 
hand side shows the case A = 2, where A = is given as reference (dashed line). Note 
the avoided crossings of the energy levels. The right-hand side shows the spectrum for 
A = 4. Here two rungs of different energy ladders calculated with the approximation 
« (dashed hues, cf. equation ijSOI)) are shown for comparison. 



These operators fulfil the commutator relations 



and 



L.N 



K . 



0. 



2k^L. 



2KL. 



(34) 



(35) 



From the commutator relations of i^^, k and L one concludes that for any function / 
of the shift operators alone, one has 

fik,k^)L = Lfi-k,-k^). (36) 

The Schrodinger equation for the time evolution operator If then reads 



4" 



6 



FdN 



U. 



In analogy to the single-band system [19] we make the ansatz 



U 



-iCit)N 



A{k,k\t) + LB{k,k\t) 



(37) 



(38) 



where C(t), A{k,k\t) and B{k,k'^,t) are yet unknown functions or operators. To 
abbreviate the notation we do not write the arguments of the operator valued functions 
A(t) and B(t) in the following explicitly. Substituting equation (|HHj) into (|H7j) leads to 

nC7(t)iVe-'^W^ (^i (t) + LB (t)) + i^e-'^^^ (t) + L^ (t) 
= FciiVe-'^(*)* (t) + LB (t)) - ^ke-'^^''>'^ (i (t) + LB (t) 
_^i^te-i^W* f i (t) + LB (t)) + h e-^^W^ f i (t) + LB {t)) . (39) 
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Here and in the following, the dot marks the derivative with respect to the time t. The 
first term of the left and the right side of equation (jH^ cancel each other if 

t 



hC(t) =Fd 



Cit) = ^ I F{r) dr. 



(40) 



The remaining terms in equation ()39p can be simplified by multiplying with e^*"'-*^^ from 
the left and using the relations 

^^cm-^^-icm =^-m)K^ (41) 
^^cmp.^^-icm ^ gicw^t ^ (42) 



which follow from the Baker-Hausdorff formula 

^g-.y ^ ^ + ^ ^ + ^ ^ ^] ] + ^ + 

Since L and N commute and LL = 1, one finally arrives at 

inX it) + ^e-'^W^i (t) + ^e^^Wi^U it) - p (t) 



(43) 



(44) 



Obviously, this equation is fulfiled if each bracket is zero. The result are two coupled 
differential equations 

A 



(45) 
(46) 



i (t) = ii| (^e-^^Wi^ + e'^Wi^t) A (t) - i-^B (t) , 
Bit) = - iA (e-^Wi^ + e'^W^t) 5 (t) - i^A (t) , 

that only depend on the commuting operators K and K'^ . With the initial condition 
f/(0) = 1 and the equations (jHHj) . (jl^j) and (jlUj) we conclude 

A 



^(0) = 1, 
5(0)=0, 



A{0) 



Ah 



5(0) = -i— 



(47) 
(48) 



Without loss of generality, we continue with the representation of k and k"^ in the 
K-basis of the single-band system 

d J2 e^"''^), (49) 




27C 



where the operators are diagonal with respect to the quasi-momentum k: 



{K\k\K) = e'"'^52. (k - k') and (/t'l^V) = e"'"'^^ 



2jT_[K — K 

d 



(50) 



with the ^-periodic delta function 62^1. The operators A(i^, i^"!", t) and B{k,k\t) 
in the K-basis are denoted by A^^it) and Bf^{t) which leads to the following system of 
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ordinary differential equations 

A. (t) = i^ cos {Kd - C{t)) (t) - i^B^ (t) , (51) 

(t) = - cos (Kd - Cit)) it) - i^A, (t) . (52) 

From now on we consider tlie case of a constant force C{t) = By substituting 

dFt d dF d 

X = Kd — and — = — (53) 

h dt h dx ^ ^ 

one obtains 

Ki^) = - i^ cos(x)A«(x) + ^^B^{x) , (54) 

B'^ix) = i^ cos(x)5^(x) + ^^A^{x) , (55) 

where the prime denotes derivation with respect to x. Differentiating again with respect 
to X decouples these equations with the result 

^"(^) + ((w)' - i (w) sinx + (2^)' cos^x) A,{x) = , (56) 
Ki^) + (( w)' + i (w) sinx + (2^)' cos^ x) S.(x) = . (57) 



Both equations are of the type of Whittaker-Hill differential equations as described, e.g., 
in [20]. A solution of the equations (|56p . (|57|) in closed form is only known in the cases 
5 = or A = 0. Needless to say that the equations and can be solved in the 
case F = 0. There always exists a solution of the differential equations (j56|l . (|57jl which 
justifies the ansatz (jHEj), even if we can not construct such a solution in closed form. 



4-2. Time evolution 

A wave packet, whose dynamics is governed by the double-periodic Hamiltonian (PJ 
is expected to show both, Bloch oscillations and Zener tunneling simultaneously. The 
tunneling rate at the edge of the Brillouin zone is approximately given by the famous 
Landau-Zener formula [21,22] 

In this chapter we study Bloch-Zener oscillations, a coherent superposition of Bloch 
oscillations and Zener tunneling. Let \'^a,n) be the eigenfunctions of the Hamiltonian 
(H)), called Wannier-Stark functions, as introduced in section 121 The index a = 0, 1 
labels the two energy ladders. We expand the Wannier-Stark functions in the K-basis, 




7r 7r 



2d 2d 

with unknown coefficients aa,n{'^) and ba^n^K). According to equation (fTTj) . the Wannier- 
Stark functions with different site indices n are related by the translation operator Tm 
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(see equation (HHI))- The Bloch waves of the two band system defined by the equations 
0, obey the relations 

T-2n|XA«) = e-^'"'^1x/3,.) • (60) 
Thus the Wannier-Stark functions ()59p can be written as 

_ZL _ZL 
2d 2d 

|^a,n)= j a„,o(«)e-''"'^1xo,.)d/€+ J b^,o{n)e-''"^''\xi,.) dn (61) 

— JL. — JL- 

2d 2d 

with projections onto the Bloch states 

(Xo,K|^a,n) = a„,o(«:) e-''"''" and (xi,.|*a,n) = &a,o(«:) e"^''^'^''. (62) 

The time evolution of an arbitrary initial state 1^/^(0)) can now be calculated in the 
Wannier-Stark basis. Expanding the initial state in the Wannier-Stark basis, 

1^(0)) = ^ Co,n|^0,n) + 5Z Ci,„|*i,„) , (63) 
n n 

the dynamics of {ip) is given by 

m)) = 5Zco,„e-^^»-*|^o,n) + J]ci,„e-^^^-*|^i,„) . (64) 

n n 

The contributions of the Bloch bands are obtained by projecting onto \Xa,K), which 
yields 

(Xo,.|^(t)) = $^Co,.e-i(2'^'^^+^»)*ao,o(«:)e-'2'^''^ 

n 

+ J2 ci,n e-^(2-<^^+'^^-^°)*ai,o(/€) e-^^"'^'^ , (65) 

n 

(Xi,.|^(t)) = 5^Co,„e-i(2"'^^+^«)*6o,o(«:)e-^^"'''^ 

n 

+ Yl e-^(^"'^^+"^-^°)*6i,o(/«) e-^2"'^'^ , (66) 

n 

where the eigenenergies flTH|) . (fT^ have been inserted and E{6, A, Fd) has been denoted 
by Eq. One observes that two Fourier series appear for the two ladders « = 0, 1, 

Ca{n + ^)=Y,c.,ne--''-<^^'^K (67) 

n 

which are vr/fi-periodic. Thus the time evolution can be written as 

(Xo,«|^(t)) = e-^^°* [ao,o(«:) Co (/^ + f ) + ai,o(«:) e'^^'^^-^^")*^! («: + f )] , (68) 

(Xi,«|^(t)) = e-^^°* [6o,o(/^) Co + f ) + &i,o(«:) e'^^'^^-^^^^^d (k + f )] • (69) 

This structure of the dynamics causes some interesting effects discussed in the following 
section. 
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4-3. Bloch-Zener oscillations and reconstruction 

The dynamics of the two band system is characterized by two periods. The functions 
Co and Ci are reconstructed at multiples of 

T.^§, ^ (TO, 

whereas the exponential function e-ji'^^~'^^o)t j^^^g period of 

27rh , , 

The period Ti is just half of the Bloch time Tb = ^ for the single-band system, 5 = 0, 
which reflects the double periodicity of the two-band model. In the case of a wave 
function confined to a single energy ladder, one of the functions Co and C\ is zero for 
all times t and the initial state is reconstructed up to a global phase after a period Ti, 
which is just an ordinary Bloch oscillation. 

In general, the functions (jH^ and reconstruct up to a global phase if the two 
periods Ti and are commensurable, 

T2 2dF r , 

The reconstruction takes place at integer multiples of the Bloch-Zener time 

Tbz = ST2 = rTi. (73) 

Examples of such reconstructing Bloch-Zener oscillations are shown in the next section. 
Furthermore we can calculate the dynamics of the occupation probability at multiples 
of Ti by the aid of the equations (jEHI), (jHTIjl . which yields 

2d 

I I (xo,. \^{nT,))\'d^ = X + Y cos vrn + , (74) 



2d 

JL. 

2d 



\{Xi,.\'^{nT,))\^dK 



, dF-2Eo 

X + Y cos — vr n + (f 

dF 



(75) 



where X and Y are real positive numbers. In the commensurable case ()72|) we find 



2d 

j |(xo,«|^(riTi))p dK = X + r cos (27r ^-n + ^ 



(76) 



2d 

TV 

2d 



J |(xi,«|^(nTi))|M«:=l- [x + Fcos(27r^n + y;)] , (77) 

_ TT 

where one recognizes the reconstruction at integer multiples of Tbz- In the case of 
incommensurable periods T2 and Ti, we obtain ip = whenever only one of the bands is 
occupied at the beginning. Otherwise one of the occupation probabilities would become 
greater than one for some n whereas the other one would become negative. 
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4-4- Oscillating and breathing modes 

The reconstructing Bloch-Zener oscillations discussed in the preceding section will be 
illustrated by some numerical calculations. The dynamics of a so-called oscillatory mode 
is shown in figure Because of the strong localization of the Wannier functions, this 
resembles closely the behaviour in real space. The initial state is chosen as a broad 



gaussian distribution in the Wannier representation, namely ~ e 



-n^/lOO 



By the choice 



of the parameters we have Eq = and we obtain a reconstruction at integer multiples 
of the Bloch time Tbz = Tb- Here and in the following, we take the Bloch time 
Tb = ^ = 2Ti of the single-band system 5 = as the reference time scale. Needless 
to say that many other reconstruction times Tbz can be adjusted in this way. The edge 
of the Brillouin zone is achieved at multiples of Ti = Tb/2, where the Bloch bands are 
close together and the probability for a band transition reaches a maximum. 





Figure 5. Dynamics of an oscillating Bloch-Zcner mode. The left-hand side shows 
IV'P versus 'space' and time. The right-hand side shows the occupation probability of 
the lower minband (dashed red line) and the upper band (solid blue line) versus time. 
The parameters are A — 80, dF — 1 and 6 ~ 6.734. 

When '^^~^p^° TT is close to 27r, the cosine in equations (f7H|) and (fTTj) changes very 
slightly after each period Ti, which illustrates the cosine-like behaviour of the occupation 
probability very clearly. This is shown in figure El for an oscillatory mode with the 
same initial condition as above. By the choice of the parameters A, dF and 5 we get 
Eq = —0.5178 and hence 

72 ^ 56 

Ti 57 ■ 

Thus we get a reconstruction at multiples of Tbz = "^STb- The slowly varying oscillation 
of the occupation between the bloch bands is similar to the Rabi oscillations for a two 
level system. 



(7^ 
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Figure 6. Occupation probability of the lower minband (dashed red line) and the 
upper band (solid blue line) for an oscillatory mode versus time. The parameters are 
A = 80, dF = 1 and 5 = 17.19. 




Figure 7. Dynamics of a breathing mode. The modulus of the wave function is 
plotted versus 'space' and time. The parameters are A = 80, = 1 and 6 = 6.734. 



Let us now consider a wave packet that is initially located at a single site (see figure 
[7j). In the case of a single-band system {6 = 0) one finds a breathing behaviour: The 
width of the wave packet oscillates strongly, while its position remains constant [12]. 
The parameters are chosen as in figure so that reconstruction happens again after 
integer multiples of Tbz = Tb- One observes that the enveloping structure with the 
period Tbz = Tb is overlayed by a breathing mode of smaller amplitude, which can be 
interpreted as a fraction of the wave packet staying in one band all the time. Therefore 
the characteristic period of this part is Ti = Tb/2. The gradient of the dispersion 
relation of a band changes at the edge of the Brillouin zone if there is no band transition. 
Therefore the direction of motion given by the group velocity changes after Ti = Tb/2 
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and the breathing amphtude of this part of the wave packet is thus smaller than the 
amplitude of the complete breathing mode. 

4-5. Beam splitters 

The Bloch-Zener transition between the two minibands can be used to construct a matter 
wave beam splitter with controllable 'transmission'. To see this, note that interband 
transitions mainly take place at the edge of the Brillouin zone (see figure Ej). In the 
case of a partial band transition, the wave packet is spatially separated, as illustrated 
in figure El 




Figure 8. Splitting of a gaussian wave packet in 'real' space. The left hand side 
shows IV'P versus 'space' and time. The right hand side shows the absolute value (red 
dashed line) and the real part (blue solid line) of the wave function at i = 0.5 Tb versus 
'space'. The parameters are A = 80, dF = 1 and i5 = 6.734. 

The time evolution of the squared modulus of a wave packet with a broad initial 
distribution ~ Q-n'^/wo jg gj-^own in figure |H1 At time t = 0.5 Tg, both parts of the wave 
packet have a gaussian profile, but the lower part of the wave function consists of states 
with opposite momentum. Hence the phase of the wave function and therefore its real 
part oscillate strongly. Usually both parts merge again after a full Bloch period Tb (cf. 
figure El). However, a permanent splitting of the wave packet can be obtained by flipping 
the sign of F at t = 0.5 Tb. Thereby both parts of the wave packet are transported in 
opposite directions (compare [23]). This effect is shown in figure El (same initial state 
as above). After flipping the field F, the external force is constant and therefore both 
splitted wave packets continue their Bloch-Zener oscillations separately. Furthermore it 
is possible to control the occupation of each branch by the choice of 6. 

The transition probability at t = 0.5 Tb for a wave packet initially located in a 
single band is in good approximation given by the Landau-Zener formula ()58|). The 
numerical results for the occupation of the different branches (compare figure IE} at 
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Figure 9. Splitting of a gaussian wave packet in 'real' space. Shown is jV'P versus 
'space' n and time t. The parameters are A = 80, d ~ 1, 5 — 6.734 and \F\ — 1. 




Figure 10. Occupation probability of the above gaussian wave packet at t ~ 0.5 Tg 
versus S. The red dashed line shows the occupation in the interval —100 < n < —40 
whereas the blue solid line shows the occupation in the interval —41 < n < 20. The 
parameters are A = 80, d = 1, and F — 1. 



t = 0.5 Tb versus 5 are shown in figure E3 Such a controlled manipulation of the band 
gap can be achieved, for example, in optical lattices. In this case it is also possible to 
construct a beam splitter in a different way, namely by 'switching off' 6 completely at 
t = 0.5 Tg. Instead of Bloch-Zener oscillations, simple Bloch oscillations will appear 
with the advantage that the two wave packets can be separated by low-loss transport 
as described in [12,23]. In any case it is possible to control the amplitude of each of the 
splitted wave packets by the choice of 6. 
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5. Conclusion and Outlook 

In conclusion, we have investigated the interplay between Bloch oscillations and Zener 
tunneling between the Bloch bands. To this end we considered a special tight-binding 
model under the influence of a static external fleld. This model satisfactorily describes 
the dynamics of a more realistic system as demonstrated in a comparison with exact 
numerical computations [23]. An additional period-doubled potential leads to a splitting 
of the ground band into two minibands in the fleld-free case. For a non-zero external 
fleld, the spectrum of the system consists of exactly two Wannier-Stark ladders. The 
dependence of the energy offset between the two ladders on the system parameters was 
discussed in detail. 

The dynamics of the system is governed by two time-scales, the Bloch period Ti 
deflned in equation (ffUjl and the period T2 of the Zener oscillation between the two 
energy ladders introduced in equation (f7T|) . If the two periods are commensurable, 
a wave packet will reconstruct periodically up to a global phase. The occupation 
of the two minibands at multiples of the Bloch period Ti oscillates sinusoidally in 
time. Furthermore it was shown that the dynamics of the system can be reduced to a 
Whittaker-Hill differential equation. Unfortunately few analytic results are available for 
the solutions of this differential equation. 

Finally some numerical examples of Bloch-Zener oscillations were given. It was 
shown that Bloch-Zener oscillations in a double-periodic potential can be used to 
construct matter wave beam splitters, e.g. for cold atoms in optical lattices. 

A further important application of a Bloch-Zener oscillation as shown in flgureElis 
the construction of a matter wave interferometer by introducing an additional external 
potential into one of the two branches of the splitted wave packet. Hence the band 
transitions act as beam splitters in a Mach-Zehnder interferometer setup. These aspects 
will be discussed in more detail in a future publication [23]. 

Appendix 

The time independent Schrodinger equation in the basis of Bloch waves consists of two 
coupled differential equations (jH)) and ((Tj). Using the ansatz 

(Xo,.|^) = Q.e-^^^/^ and = P.e^''^/^ (79) 

we will show that in the case of a two-band system there exist exactly two energy 
ladders. A similar proof can be found in [24] for a different system. Substituting (fTUjl 
into equations (jHl) and ((Zj) we obtain 

d_( Q.\^( ^o,« - dF M: 

The dispersion relations -Eo,k' reduced transition matrix element Mf^ are 

7r/(i-periodic. Therefore we can apply Floquet's theorem, which says that a system of 
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linear homogeneous differential equations with periodic coefficients as given above is 
solved by a fundamental system of the structure 

X, = Y,e'^^, (81) 

where is a Tr/rf-periodic 2x2 - matrix and Z is a constant complex 2x2 - matrix. 
Let the eigenvalues of Z, called characteristic exponents, be denoted by Zq and Zi. They 
are only defined up to multiples of 2d. Multiplication of the fundamental system (jHTjl 
with the eigenvectors of Z from the right leads to 




QK,n 
P 

with 7r/(i-periodic functions /k,„ and (^^.n- Here, n = 0,1 distinguishes both possible 
eigenvectors. Because (xo,k|^) and (xi,k|^) have to be vr/ci-periodic, too, the phases 
are given by 

TT E IT 

- 1-— + i-Zn = -il-Km (83) 

with m ^J_. The energies are therefore 

En,m = 2mdF + FZn . (84) 

Obviously, there exist exactly two different energy ladders as long as the eigenvalues of 
Z are distinct. Due to the fact that the Hamiltonian is hermitian, the energies -En,m , 
and thus the z^, are real numbers. 
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